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Abstract 

The results of the computer hunt for the primes of the form q = m 2 + 1 
up to 10 20 are reported. The number of sign changes of the difference 
Ttnix) — % fn r-^ u i ^ and the error term for this difference is investi- 
gated. The analogs of the Brun's constant and the Skewes number are 
calculated. An analog of the B conjecture of Hardy-Littlewood is for- 
mulated. It is argued that there is no Chebyshev bias for primes of the 
form q = m 2 + 1. All encountered integrals we were able to express by 
the logarithmic integral. 
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1 Introduction 

This paper is devoted to investigation of the set of prime numbers 

Q = {2, 5, 17, 37, 101, 197, 257, 401, 577 . . .} (1) 

given by the quadratic polynomial m 2 + 1 and let q n denote the n-th prime of this 
form. By the conjecture E of Hardy and Littlewood [19] the numbeiQ7r g (a;) of primes 
q < x of the form q = m 2 + 1 is given by 



where 

C * = 11 ( 1 ~ i — ) = 1-372813462818246009112192696727 ... (3) 

p>3 VP/ 

The primes q n were investigated in the past both theoretically and numerically. One 
of the strongest theoretical results is the theorem of H. Iwaniec [22] , who proved that 
there exist infinitely many integers m 2 + 1 which are 2-almost-primes. In 1998 H. 



1 we adopt here the convention that all functions describing usual prime numbers will have the 
subscript q when applied to primes of the form q — m 2 + 1 



2 



1 Introduction 



Iwaniec and J. Friedlander [T5] have proved that there is an infinity of primes given 
by the polynomial of two variables m,n of the form m 2 + n 4 , thus the case of the 
polynomial of one variable m? + 1 is not covered by their theorem. Quite recently 
there appeared two papers by S. Baier and L. Zhao [21 E], treating the more general 
problem of the primes of the form m? + k in average over square-free parameter 
k in appropriate intervals. In the computational part we should cite the papers 
by Shanks [3B1 EE] and Wunderlich [32]; in the last paper the table of ir g (x) for 
x < 1.96 x 10 14 is given. 

By analogy with the case of all primes, where substitution of the logarithmic 
integral Li(x) = duj log(w) instead of x/log(x) gives better approximation for 
7r(x), we recast the original Hardy and Littlewood conjecture E in the form: 

n q (x) ~ -C q f ^ (4) 
qK ' 2 9 J 2 v^logn 1 ; 



/ y/x 2jx 8^x 48^ 2 n ~\n - l)\y/x 

~ 1 I" i 2 ^ i 3 ^ i 4 H • • ■ H : n h ■ 

\logx log x log x log x log X 

The series in parenthesis above is asymptotic one: the terms initially decrease, 
but for sufficiently large n they become to increase. Thus this series has to be cut at 
such no, which depends on x, that the nn-th term (and consecutive terms) is larger 
than previous one with n (x) — 1 what gives for the threshold at which the series 
Q should be cut the inequality: 

n Q (x) > -log(ac) + 1. (5) 

Besides this asymptotic representation the integral Q can be linked to the logarith- 
mic integral by the change of variables u = t 2 : 



= /^ = li(V&)-li(>/S). (6) 

J a V^\0gU jr~ a log* 



-Vb 

Here we use the following convention for the lower limit of integration: 







Integration by parts gives the asymptotic expansion: 

,., . x x 2x 6x n\x 

~ TR^X + RRR"R + RRR"R + RRRR + ' ' ' r~nTO' ( 8 ) 



log(x) log 2 (x) log 3 (x) log 4 (x) log 



which should be cut at n Q = log(a;). There is a series giving li(x) for all x and 
quickly convergent which has n\ in denominator and log n (x) in nominator instead 
of opposite order in ([8]) (see [U Sect. 5.1]) 

log^x) 
n ■ n\ 



^ lOE (X) 

li(x) = 7 + log log(x) + , for x > 1 , (9) 



n=l 
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Here 7 = 0.57721566490153286... is the Euler-Mascheroni constant. Even faster 
converging series was discovered by Ramanujan 0, p. 123]: 



=7 + log(log(x)) + v ^^ ( 1)W "j 2 ( y )f £ far x> 1.(10) 

n=l ' fc=0 

Skipping the number 2 = l 2 + 1 all odd primes q n can be expressed by the form 
m 2 + 1 only for m even, thus we put m = 2k, i.e. we are looking for primes of the 
form 4k 2 + 1. We have collected in one file of the size roughly 3.5 GB values of 
k (thus the prime 2 = l 2 + 1 is absent here!) up to k = 4999999978 what corre- 
sponds to all primes of the form m 2 + 1 < 10 20 . The compressed data occupies 760 
MB and is available for downloading from j http : //www, if t .uni . wroc .pl/~mwolf /| 



4k-2+l-data.zip, Initially up to 6.65 x 10 16 I have used my own program written 
in Fortran for Alpha DEC workstation, but to reach 10 20 we have used the free pack- 
age PARI/GP [39J developed especially for number theoretical purposes. To scan 
the interval (6.65 x 10 16 , 10 20 ) it took about 10 days of CPU time on the PC with 
the clock 2.5 GHz using the built-in PARI very fast function isprime(p) . The Table 
I gives comparison of ir g (x) with conjectures ^ and Q. Among those 312,357,934 
values of k there were 11,864,645 such k that they in turn were primes, i.e. there 
were 11,864,645 pairs of numbers (p, 4p 2 + 1) both being prime. We have checked 
separately that up to 10 22 there are 96,817,209 such pairs (p, 4p 2 + 1) [40] . 

As it is seen from the fourth column, the actual number ir q (x) is always larger 
than prediction (J5]). Contrary to this the ratio of n q (x) to the integral Q sometimes 
is larger than 1 and sometimes is smaller than 1, see the last column. It means, 
that the difference n q (x) — \C q j* ^" gM changes the sign, see the Sect. 3. In the 
next Section 2 we discuss the problem of the error term, in Sect. 4 the analog of 
the Brun's constant is calculated. Sections 5 and 6 contains some heuristics about 
the distribution of k in 4k 2 + 1 giving the prime. We formulate an analog of the B 
conjecture of Hardy and Littleewood for the case of primes q = m 2 + 1. In Sect. 
6 we discuss the distribution of the gaps between consecutive k's giving the prime 
4/c 2 + 1. Heuristic arguments allow us to make conjecture about the growth of the 
difference q n+ \ — q n = (9( v /g^log 2 (g n )). The last Section 7 contains the discussion 
of the analog of the Chebyshevs Bias for the case of primes in Q. 



2 The problem of error term 

Nothing is known about the error term for the formula Q (see however [21 E]), thus 
the only way to gain some information and intuition is to appeal to the available 
computer data. The Figure 1 presents the plot of the difference 



|A,(x) 



7T 9 (x) - ~Cg 



du 



ulogu 



for x G (10, 10 20 ). As it is seen from this plot the graph of the error term is very 
erratic, thus the plot of the maximal value of the absolute difference: 



u(x) = max |A„(t)|, 

2<t<x H 



(12) 
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2 The problem of error term 



what is a kind of envelope for |A g (x)|, is also plotted in green in the Figure 1. 

The error term present in the ordinary Prime Number Theorem (PNT) under 
the Riemann Hypothesis is ^/xlog(x)^. 

ir(x) =U(x) + 0(y/xlog(x)). (13) 

which can be written in the slightly weaker form 

tt(x) = Li(a;) + 0{x 1 i +e ). (14) 



TABLE I 



X 


7T q (x) 


Cqs/xj log(x) 


Ti q {x)/ eq. (2) 


formula (4) 


ir q {x)/eq. (4) 


10 6 


112 


99 


1.12713 


122 


0.91869 


10 7 


316 


269 


1.17325 


318 


0.99440 


10 8 


841 


745 


1.12847 


855 


0.98321 


10 9 


2378 


2095 


1.13516 


2357 


1.00888 


10 10 


6656 


5962 


1.11639 


6610 


1.00696 


10 11 


18822 


17140 


1.09815 


18787 


1.00184 


10 12 


54110 


49684 


1.08909 


53971 


1.00258 


10 13 


156081 


145028 


1.07621 


156386 


0.99805 


10 14 


456362 


425861 


1.07162 


456405 


0.99991 


10 15 


1339875 


1256912 


1.06601 


1340089 


0.99984 


10 16 


3954181 


3726285 


1.06116 


3955222 


0.99974 


10 17 


11726896 


11090399 


1.05739 


11726340 


1.00005 


10 18 


34900213 


33122538 


1.05367 


34903278 


0.99991 


10 19 


104248948 


99229889 


1.05058 


104251624 


0.99997 


10 20 


312357934 


298102838 


1.04782 


312353427 


1.00001 



The ^Jx behavior is confirmed by the computer data, see e.g. [32} Table 14, p. 
175] or [17, Table 5 and 6], where the difference Li(x) — -nix) has roughly half digits 
of the value k{x). Because there are roughly ^Jx candidates for primes of the form 
m 2 + 1 up to x it is natural to expect that the error term for Q will be square root 
of the error term for PNT. This heuristic seems to be confirmed by the fact, that 
u(x) is well approximated by the power-like error term: 

aix h , ai=0.38..., fa = 0.22 . . . . (15) 

and indeed here {3± fa 1/4. This function was obtained by fitting the straight line to 
the points log(o;(x)) vs log(x) for x > 10 9 by the least-square method and to bound 



the difference u(x) from above it is sufficient to shift the above curve (15) parallel 
up to leave the plot of |A ? (x)| below. In the Figure 1 the function 5xf was chosen, 
but at least for x < 10 20 the smaller choice for the constant hidden in the big-O in 
u(x) = 0{x 131 ) will also do. These heuristic arguments and computer data lead us 
to guess the following 
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Conjecture 1: 



du 



ulogu 



+ 0(x 



(16) 



More stringent error term 0(x* v/log(x)) is also a possibility which cannot be ruled 
out by available data. Let us notice, that (16) is heuristically supported by the 
relation ([6j and ir(y/x) = \i(^fx) + 0(x 1 ^). Another support in favor of (16) will be 
given in Sect. 7 (see Fig. 8). 

The difference A q (x) fluctuates roughly symmetrically around zero. As the com- 
puter check of possible future oscillations theorems we present in the Fig. 2 the plot 
of A q (x) / x 1 ^ . The amplitude of this function practically does not change in the 
interval (10 2 , 10 20 ) giving support in favor of the error term 0{x 1 ' i ). 



3 An analog of the Skewes number for primes of the form 

m 2 + 1 

It turns out that the there is a lot of sign changes of the difference 

A q (x) = n q (x) - \c q f f\ (17) 

in the investigated interval x G (2, 10 20 ). In the generic problem of all prime numbers 
it was shown by J.E. Littlewood in the 1914 [28] (see also [13]) that the difference 
between the number of primes smaller than x and the logarithmic integral li(x) 
infinitely often changes the sign. The smallest value xs such that for the first time 
the difference A(x) = ir(x) — li(x) changes the sign is called Skewes number and 
the lowest present day known estimate of the Skewes number is around 10 316 , see 
[5] and [UJ. However in the case of the primes given by the quadratic polynomial 
m 2 + 1 the first sign change of the difference A q (x) occurs already at the prime 
qi3 = 2917 = 54 2 + 1 and there are 20634 such sign changes up to 10 20 . Let v q (T) 
denotes the number of sign changes of the function A q (x) for x G (2, T). The Fig.3 
presents the plot of the function v q (T) for T < 10 20 . The fitting of the power- 
like dependence of v q (T) on T gives parameters which depend on the number of 
discarded initial points. For example fitting \og{y q (T)) vs log(T) for T G (10 7 , 10 20 ) 
gives u q (T) ~ T - 23308 while for T G (10 12 , 10 20 ) we obtained u q (T) ~ T 23834 , and 
we have checked for other intervals of T that the first digits 0.23 persists, thus we 
write: 

u q (T) w a 2 T^ 2 a 2 = 0.38 . . . , fa = 0.23 . . . (18) 

Let us mention that for the case of all primes Knapowski [25] proved that the number 
of sign changes of A(x) = in the interval (1, T) 

v(T) > e~ 35 log log log log T (19) 

provided T > expexpexpexp(35). There is a remarkable coincidence in the values 



of the parameters a and (3 present in the fits (15) and (18) and when the functions 
uj(x) and v q (T) are plotted on the same graph they appear very close to each other, 
despite the fact that u(x) and v q (T) represent quantities at first sight unrelated. 
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4 The analog of the Brun constant 



4 The analog of the Brun constant 

Because Y^m=i 7? = 71-2 A> < oo thus the sum of reciprocals of all primes of the form 
q = m 2 + 1 is trivially convergent: 

1 1 1 1 1 1 

> - = - + - + — + — + + ...<oo (20) 

^ q 2 5 17 37 101 K J 

but the actual numerical value of this sum is unknown [31 J . In 1919 Brun [9] has 
shown that the sum of reciprocals of all twin primes is finite: 

*-G+SMI+Hii+n)+-<- (21) 

Numerically B 2 = 1.9021605823 . . ., see [2"9"j . [8]. It is natural to call the above sum 
(20) the Brun's constant for primes of the form q = m 2 + 1 and denote it by B q . 
From the computer data we can calculate the finite size approximations: 

B q (x) = £ I. (22) 

x>qeQ 

From the integral test for convergence of the series in the form: 



oo 

E 

n=N 



f(n) < f(N) + / /(«) du. (23) 



N 



we have: 



B q 



El 1 ^ , s sr^ 1 ^ , s 1 f°° du „ . . 1 1 

-+52 -<B q {x)+52~ 2 <B q {x) + - + / ^=B q (x) + - + ^ 

(24) 

Thus from this trivial inequality we can expect for x = 10 20 the accuracy of 10 digits 
for B q . Indeed, from the second column of the Table II we see that the number of 
stabilizing digits of B q (x) is roughly half of the digits of the exponent of x in the 
first column. However, using some heuristics it is possible to obtain from i3 9 (10 20 ) 
15 digits of B q . Namely, we can obtain the analytical formula for dependence of 

it follows that the chance to find a prime of the 
thus we can write: 



B q (x) on x. From the equation (|4 
form m 2 + 1 around x is 



2y / Slog(o;) 



x>q£Q 

Integrating by parts gives: 



w = E q * - f i 



du 



u 3 / 2 log(u) 



(25) 



u 



2 + * + + (26) 



3 / 2 iOg(w) y/x\ogX A/xlog 2 ^ y/x\og n X 



2 Let us remark, that in the case of all primes tt(x) ~ J 2 lo g" M ) = \ g( x ) + ■ • ■ a s well as of twin 
primes 7T2(x) ~ C2 f% i Q ^" u ) = C2 ] og ? ^ + • ■ • the chance to find a prime around a; following from 
dividing the first terms on r.h.s. by x coincide with integrands on l.h.s, what is not true in the 
case of |4|), where the factor | makes the difference. 
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This series is asymptotic one and the condition for the dropped terms is n > \ log(x) 
- the same as threshold ([5J . But by the change of the variable u = 1/ s/i it is 
possible to express the above integral by the logarithmic integral: 



du 



a W 3/2 log(w) 



1 

7h 



-ii 



(27) 



This formula is useless for our purposes, because ^ or (10) is valid for x > 1. The 
analog of (25) for usual Brun's constant is given by 



where C> 



2C 

B 2 (oo) = B 2 (x) + — 2 - 
log(x) 

) = 0.66016 ... is the Twins constant. 



(2f 



p>2\ J - (p-1) 2 ' 

The third column in Table II gives the sample of values 



B*(x) 



B q {x) + ^ 



du 



m 3 / 2 log(u) 



(29) 



which are supposed to be constant and equal to B g (oo). As it seen from the Table 



II indeed with increasing x growing number of digits of the sum (29) is stabilizing. 
To produce the data for this Table we calculated in PARI the value of C q with 
over 30 digits accuracy using the formula (10) from |36j. We calculated the finite 
approximations B q (x) in DEC Fortran using the quadruple precision (REAL* 16) 



with 33 decimal digits and the Mathematica v. 7 was used to calculate integrals (26) 
with over 30 digits of accuracy. 

TABLE II 



X 


B q (x) 


B*(x) 


10 2 


0.79575154653780163 


0.81798411614326626 


10 3 


0.81119372199008314 


0.81625948818529443 


10 4 


0.81335296609757082 


0.81460891435244917 


10 5 


0.81432340308206016 


0.81465046053349854 


10 6 


0.81450766696668914 


0.81459563341678427 


10 7 


0.81457232824055968 


0.81459653781164721 


10 8 


0.81458971836435488 


0.81459649673657666 








10 13 


0.81459655805079256 


0.81459657169340710 


10 14 


0.81459656768254395 


0.81459657170479096 


10 15 


0.81459657051185242 


0.81459657170321129 


10 16 


0.81459657134862317 


0.81459657170292293 


10 17 


0.81459657159724819 


0.81459657170299049 


10 18 


0.81459657167131276 


0.81459657170297237 


10 19 


0.81459657169347024 


0.81459657170297623 


10 20 


0.81459657170012661 


0.81459657170298816 



It is seen from the last column that starting with x = 10 16 all first 14 digits 
remain the same — the change appears at the 14-th place after the dot. In the 
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5 Analog of the conjecture B of Hardy-Littlewood 



paper [3U] Nicely has performed complicated statistical analysis to get the 95 % 
confidence interval for the value of 25 2 . In our case it is possible to estimate the 



error appearing in (29) by using the form of the function u(x) given by Conjecture 
1. Namely, the "density" of the error for the chance 2y /E\og(x) ^° nnc ^ ^ ne P rmie °f 
the form m 2 + 1 around x is less than 0(x~ 3 ^), thus we have: 

\B q (oc)-B;(x)\ = 0(x- 3 ^). (30) 

From this we see, that for x = 10 20 the value of B q (oo) lies in the interval of 
approximate length 10~ 15 around i3*(10 20 ) and we can claim that with 15 digits 
accuracy 

B q (oo) = B q = 0.81459657170299 .... (31) 



S. Plouffe has checked using his Symbolic Inverse Calculator (http:/ /pi.lacim.uqam.ca/eng/), 
that this constant can not be expressed by other mathematical constants [31], thus 
the value of B q could be treated as a new mathematical constant. The compari- 
son of numbers in the second and third column reveals that addition of the term 
% f°° i/?^ i \ causes that about 3-4 digits more than in the values of B a (x) alone 
settle down, thus the rate of convergence of B*(x) is a few orders faster than that 
oiB q (x). 



The correctness of the choice of | in front of the integral in (25) can be checked 
by comparing the values of the equation 

C f X2 dn 

following from (25) with the actual computer data. For example for x\ = 10 6 and 
Xi = 10 20 for the l.h.s. of (32) from the computer data we get 0.000177. . . while the 



r.h.s. is equal to 0.000176 . . .. 

Let us mention that two first primes from Q give contribution | + | = 0.7 to B q , 
i.e. 86% of the total value 0.8145966 . . .! 

We can define the Brun's measure of the set of numbers S = {ai, . . . , } as the 
sum 

M *W = E rV ^ 

Mi 

i 

provided it is finite. Thus we can say that the Brun's measure of the set of twin 
primes is 1.9021605823..., while the Brun's measure of the set Q is M.b(Q) — 
0.8145965717.... 

5 Analog of the conjecture B of Hardy-Littlewood 

The conjecture B of Hardy-Littlewood [12] says that the number 

7i(d; x) = 1 =| {p < x : p and p — d are primes} \ (34) 

p<x 
p—p'=d 

of prime pairs p — d,p < x separated by d (d — 2, 4, 6, . . .) is given by 

^>^n(^)f^, (35, 

P \d 
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where 



C 2 = 2 TT ( 1 - - — ^— - J = 1.320323631693739 . . . (36) 

p> 2 V (p-i)v 



is the Twin constant. The integral in (35) again can be expressed by the logarithmic 



integral and thus calculated quickly from the series (\9\j or (10): 



b dt f b dt a 6 _ ,,s , . a b 

a io?w ~ j a ^m + M^y " M6) " lU ~ 1(Gj + bgw " M&) ( j 

The values of the gaps between primes m 2 + 1 and (m — d) 2 + 1 grow linearly 
with m, but we can formulate an analog of the conjecture B of Hardy-Littlewood 
when we will focus on the gaps between k's appearing in Ak 2 + 1. Thus let us define: 



7i q (d;x) =| {2k < yfx : 4k 2 + 1 and A{k — d) 2 + 1 are both primes} | (38) 

In contrast to all primes here d = 1, 2, 3, . . .. The Fig. 4 presents the plots of Tr q (d; x) 
obtained from our computer data for x = 10 10 , x = 10 12 , . . .x = 10 20 and d < 300. 

There is a heuristic procedure of Bateman and Horn |4] (see also [331 chap. 3]) 
allowing to guess the formula for the number of constellations of primes of differ- 
ent types. Let f = {fi(x), fz(x), . . . , fi(x)} be the set of distinctive irreducible 
polynomials with integral coefficients and positive leading coefficient such, that 
f(x) = fi(x)f2(x) . . . fi(x) has no fixed divisor > 1. Let 7r(f; x) denote the number 
of positive integers n < x such, that all fi(n), f2 (n),..., fi(n) are simultaneously 
primes. Then the Bateman-Horn conjecture reads: 

7r(f;x) -TTt K7 / —, — > (39) 

where w{p) is the number of distinct solutions to fx(u)f2(u) ■ ■ ■ fk(u) = (mod p) 
with u G {0, 1, . . . ,p — 1}. In our case fi(u) = Au 2 + 1, f 2 {u) = A(u — d) 2 + 1. It is 
well known [20, chap. VII] that equation 4u 2 + 1 = (mod p) has two solutions for 
primes of the form p = 4n + 1 and does not have solutions for primes of the form 
p = An + 3. For the case p = An + 1 if p \ d then equations fx(u) = (mod p) and 
/ 2 (w) = (mod p) have the same solutions, thus w{p) = 2. When p \ d then there 
are two possibilities: w(p) = 3 when fi(u) = (mod p) and f2(u) = (mod p) 
have one common solution, or w(p) =4 when fi(u) = (mod p) and f2{u) = 
(mod p) have distinct solutions. The equations fi(u) = (mod p) and f2(u) = 
(mod p) can have one common solution only when p \ ±2u* + d, where u* is the 
solutions of 

Au 2 + 1 = (modp). 
The solutions of this equation are of the form (see [20], p. 88]) 

v — 1 

2ux 2 = ±u*, u* ' 



Because the conditions p \ ±2u* + d for the case w(p) =3 can be written as the one 
condition p \ d 2 — Au* and there is an identity 4m* 2 + 1 = (mod p) we can write 
the condition for w(p) =3 simply as p \ d 2 + 1. Thus we have for p = An + 1 
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5 Analog of the conjecture B of Hardy-Littlewood 



[4 if p \ d, p\d 2 + l 
w{p) = I 3 if p\d 2 + l (40) 
[2 if p I d 

For p = 2 and for p = 4n + 3 there are no solutions, i.e. w(p) = 0, hence finally 
the product appearing in the Bateman-Horn conjecture takes the form: 



w(p)/p 



1/p) 



2 



(41) 



4 TT 1 TT 1-4/p TT l ~ 2,p n 1-3/p 

J-J- (i-1/p) 2 11 (i-1/p) 2 11 (i-1/p) 2 11 (i-1/p) 2 ' 

p=3 (mod 4) V ' ^' p=l (mod 4) V ' ^' p=l (mod 4) V ' ^' p=l (mod 4) V 1 ^' 

We can get rid of the product over p j d by extending it to the product over all 
p = 1 (mod 4) and simultaneously by dividing by an appropriate term. Because the 
conditions p \ d and p | d 2 + 1 cannot be satisfied simultaneously these additional 
factors can be incorporated into the last two products above. Finally the factor 
describing oscillations takes the form (we separated 4 to cancel it later with 4 coming 
from the degrees of fi(u) and /2(w)): 



W{P)/P -4C lP (d), (42) 



(1 - l/pY 
where the constant C\\ 

Cl= II 7^~W II P (^~rv = 0-975245556223143537223292783 . . . 

(43) 



(p-i)2 ii (p-iy 

p=3 (mod 4) v ' p=l (mod 4) v ' 



and -P(d) denotes the product: 

*w>= n j£ n m 

p=l (mod 4) P=l (mod 4) 

p|d p|rf 2 +l 

In above expressions the condition p = 1 (mod 4) means that products are over 
primes p > 5, thus all these products are positive. In (44) the conditions p \ d and 
p | d 2 + 1 are fulfilled only by finite number of p's, hence it is obvious that these 
products are convergent. 

Finally we obtain the number of such k < x that both fi(k) = (2k) 2 + 1 and 
f2{k) = (2(k — d)) 2 + 1 are prime: 

»</.,*;*) =c n n —J'v^ < 45 ' 

p — 4 - LJ - p — 4 /i log (2m) 

P e1 (mod 4) r P e1 (mod 4) r ^ 1 6 V > 

P\d p\d 2 +l 



We put here for a while 1 as the lower limit of integration, since k — 1 gives the 
prime 4 ■ l 2 + 1 = 5 (let us notice, that Ramanujan often did not specify the lower 
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limit of integration, see [51 p. 123]). Because we skip 1 in m 2 + 1 in manipulations 
of integrals below, alternatively we can say that the lower limit of integration is 
V / 5/4 = 1.118033989.... 

Usually we are interested directly in the number of primes Ak 2 + 1 < x and after 
the change of the integration variable u = y/i/2 we have 



^ 2 du r b dt 

log 2 (2«) J a Vi\og 2 (t) 



(46) 



and finally for the quantity 7r q (d;x) defined in (38) we obtain Conjecture 2: 



/ i \ ^ TT P~ 2 TT P~ 3 f dt 

TT a {d\x) = C\ / — = k h error term (47) 

11 p-4 11 Jt\og 2 (t) 

P e1 (mod4) r p=l (mod 4) ^ J 5 V blu &\ b ) 

p\d p\d 2 +l 

We cheated a little here replacing 4 by 5 as the lower limit of integration; better pos- 
sibility as the lower limit of integration is perhaps the Soldner-Ramanujan constant 
fj, = 1.45136923488338105... defined by li(/i) = 0, see p.123]. Appearing here 
integral by the change of the variable t = u 2 can be expressed by the logarithmic 
integral: 

/ r d \ = -(H^b) ~ li(v^)) + r-^r " r-^r- (48) 
J a Vi\og 2 (t) 2 K K } KV ,J log(a) log(6 1 ; 



The product P(d) (44) is responsible for characteristic oscillations seen in the 
Fig. 4. The conjecture 2 agrees with the computer data quite well. Instead of 
producing some table to corroborate this statement, we give "visual argument": 
in the Fig.4 by the red boxes are plotted values of the quotient ir q (d;x)/P(d) for 
x = 10 12 ,...x = 10 20 and d < 300. In green are plotted values of the integral 
f* du/ y/t log 2 {t) for the same values of x. 

6 Heuristics on the gaps between adjacent k 

It is interesting to restrict the analysis from the previous section to the case of 
consecutive values of k giving the prime 4A; 2 + 1. Let us define the quantity: 

h(d; x) = {number of pairs k < k! = k + d, such that Ak 2 + 1 and Ak' 2 + 1 < x are 
consecutive primes of the form m 2 + 1} = 1. (49) 

qn-qn-i=4d(y/q n -l-d) 

The Fig. 5 presents the plot of h(d; x) obtained from our computer data for 
x = 10 10 ,x = 10 12 ,...x = 10 20 . On the semi-logarithmic scale the points display 
characteristic oscillations around straight lines representing the fits to exponential 
decrease obtained by the least-square method. To describe the oscillations we take 



the product P(d) from (44). The Fig.5 presents the corroboration of this mecha- 
nism of oscillations: dividing the values of h(d, x) by the product P(d) leaves pure 
exponential decrease — small deviations could be attributed to the fluctuations of 
h(d, x) and incorporated into the (unknown) error term. 
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6 Heuristics on the gaps between adjacent k 



The Figure 5 suggests the following 
Ansatz: 

h(d, x) = C 1 P(d)B(x)e- dA(x \ (50) 

The functions A(x) and B(x), giving the slopes and the intercepts of straight lines 
seen in the Fig. 5, can be determined by exploiting two selfconsistency conditions 
that h(d, x) has to obey just from the definition. First of all, the number of all gaps 
between k's is by one smaller than the number of primes of the form 4k 2 + 1 smaller 
than x: 

K(x) 

h(d,x) = 7T q (x) — 1, (51) 

d=l 

where K(x) denotes the largest gap between two consecutive k, k! < x. The second 
selfconsistency condition comes from the observation, that the sum of distances 
between adjacent k is equal to the k producing the largest prime q = 4k 2 + 1 < x. 
For large x we can write: 

K(x) 

Y,h(d,x)d=Y- (52) 

d=l 

The erratic behavior of the product P(d) is an obstacle in calculation of the above 



sums (51) and (52). Thus we will replace P(d) by the mean value: 

n 

J2 P ( d ) =ns + E{n), 



(53) 



d=l 



where we assume that the unknown error term E(n) is an increasing function of n 
which grows slower than n: 

li m = o (54) 



n— >oo n 



what means that: 



1 n 
lim — y 



n 



p- 2 



n 



p — 3 



p — 4 J-- 1 - p — 4 

p=l (mod 4) P=l (mod 4) 

\ P\d p|rf 2 +l ' / 



(55) 



There is known at least one example of the error term which grows slower than n for 
the similar problem. Namely, E. Bombieri and H. Davenport [7] have proved that 
the number 1/ ri P >2(l — (p-i) 2 ) ^ s ^ e ar it nme ti ca l average for the product Yl P \d ^E^'- 



E n 

d=l p\d,p>2 



p — 1 



71 



+ o(io g 2 H). 



(56) 



p>2^ 



(P-I) 2 



To get rid of P(d) in (51) and (52) the Abel summation formula can be used in 
the form: 



n— 1 



^2 a i b i = ~ y]S(i)cj + A(n)b r: 



i=l 



i=l 
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where S(i) — ai + . . . a, and q = b i+ \ — q. Putting here Oj = P(i), hi 
replacing E(l) < E(2) < . . . < E(n — 1) by larger E(n) we obtain: 



n-1 



Y,P(l)f(l) = sJ2f( l ) + 0(f(n)E(n)). 



f{i) and 



(57) 



i=i 



In our case /(/) = C\B{x)e A ( x > 1 for equation (51) and /(/) = C\lB(x)e A ^ x ' 1 for 
equation J52J) . From the Fig.5 we see, that h(d, x) decreases exponentially with d and 



to solve fl51 



and (52) for e A ^ and B{x) we have to drop the term 0(f(n)E(n)). 



The sums (51) and (52) are the geometrical and differentiated geometrical series 



respectively; because h(d, x) decreases exponentially with d we have replaced K(x) 



in (jSTj) and ( |52| ) by oo and (plj) , (|52|) turn into the equations: 

7T q (x), 



sC 1 B(x)e^ A< -^ 



1 _ e -A(x) 

sC l B(x)e- A{x ^ 



x 



(1 - e- A W) 2 2 ' 
The solutions for e~ A ^ and B(x) of the above equations are: 

27T n (x) 



(58) 
(59) 



-A{x) 



1 - 



X 



-, m 



B(x) 



2<(x) 



8 CiVx(l - ^ 



(60) 



(61) 



Hence we have finally: 



h(d; x) 



Sy/x 



27T n (x. 



X 



d-1 



_l_ error term. 



(62) 



The Table III gives a comparison of the formulas (60) and (61) with the slopes 



and intercepts of \og(h(d; x)/CiP(d)) vs d obtained from the computer data by 
means of the least square method (roughly 1/4 values of h(d;x) for largest d were 
discarded to avoid the fluctuations of h(d; x) in the region of large d). 



TABLE III 



X 


e -A{x) 


I 2tt,(x) 


2ir q (x) 

e v^s 


B{x) 


1 2tt2(x) 


2tt2(x) 




sCi ^/x—2ir q (x) 


sCly^X 


10 15 


0.91014186 


0.91525886 


0.91875009 


75421.77 


65826.02 


60247.85 


10 16 


0.91611160 


0.92091638 


0.92396266 


206903.63 


180179.71 


165930.45 


10 17 


0.92141730 


0.92583260 


0.92851624 


570171.11 


498478.78 


461507.91 


10 18 


0.92632759 


0.93019957 


0.93257992 


1573509.35 


1389610.28 


1292614.89 


10 19 


0.93060739 


0.93406718 


0.93619375 


4388166.12 


3904615.30 


3647172.98 


10 20 


0.93444100 


0.93752841 


0.93943975 


12318241.09 


11044184.29 


10354236.57 
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6 Heuristics on the gaps between adjacent k 



The expression for the mean value s of the product P(d) can be obtained in the 
following way: Because the pairs of primes of the form (2k) 2 + 1 and (2k + 2) 2 + 1 
(Shanks calls in [37] such pairs Gaussian Twins) correspond to d — 1 thus they 
are necessarily consecutive (q n , q n +i) and the formula for the number h(l, x) of such 



pairs smaller than x obtained from (62) has to be equal to 7T g (l; x) from (47): 



h(l; x) 



2*j(s) 
Sy/x 



fit 



v^log 2 (t) 



(63) 



For large x we have 7r 9 (x) = C q \fxj log(a;) and 

dt 1 



x 



l v/tlog 2 (t) 2 h( ^ } log(x) 



2^c 
\og 2 (x) 



+ ... 



where we have used two first terms of the asymptotic expansion (|8l and fortunately 



the term yfx/\og(x) cancels out leaving on both sides of (63) the same dependence 
on x and thus we obtain s 



C 2 JC X : 



lim 



n 



( 



\ 



P 



TT - 

p — 4 A A p 

p=l (mod 4) P=l (mod 4) 

p\d p\d 2 +l 



n 



(64) 



In the full (mysterious) form it reads: 



lim 

n— >oo n 



1 n 

-Y, 



( 



\ 



n 



p 



n 



p 



p — 4 J-- 1 - p — 4 

■ p=l (mod 4) P=l (mod 4) 

\ p|d p\d 2 +l 



n 



p>3 



(_l)(P~l)/2 

P 3 ! 



n 



P(P~4) 

p=3 (mod 4) (p-1) 2 llp=l (mod 4) (p-1) 2 



(65) 

We were not able to prove this identity analytically — usual methods of calcu- 
lating the sums of arithmetic functions, see e.g. [231 Chap.l], are not applicable 



here because P(d) is not a multiplicative function. The computer checking of (65) 



for large number of terms on the l.h.s. is also difficult because the calculation of 



the average has almost cubic complexity in n (i.e. finding the value of l.h.s. of (65) 
involves 0(n 3 / log 2 (n)) operations). We have calculated the sum of P(d) for d up 
to 150000 and we obtained P ( rf )/ 150000 = 1-93242674. . ., while the r.h.s. 



of (65) is (1.3728134628) 2 /0.97524552 = 1.93245368, thus the first 5 digits are the 



same. 



In [37], [HJ P- 90], [35] heuristically the formula for h(l;x) was obtained in the 
form: 



h(l,x) = F 



log (x) 



_l_ error term 



(66) 



where: 



7T 



n 



p=l (mod 4) 



P+l 
P J \p — 1 



1.9504911124462870744465855658 .... 

(67) 
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In [33] the 50 digits of this constant are given. Therefore we have C\ = F/2, s 
2Cg/F and the combination on the r.h.s of (65) can be transformed to the form: 



2C£ = 1 tt (P-2) 2 ( P -1) 2 n fp + l\ 
F 4 H p(p-A)(p+l) 2 11 U-U 

p=l (mod 4) ' p=3(mod4) VK 



Finally we state the 
Conjecture 3: 



FPi^ixlf 2n q (x)\ d - 1 
h{d, x) = — — I 1 j=— ] + error term. (69) 



For large x we can simplify considerably the above formulas by writing e A ^ 
have: 



27Tg£r) _ ^ and -B (x) = 2 ^ x )_ therefore in the limit of large x we 



h(d,x) = ir g (x)e ^ + error term. (70) 

C q y/X 

The Table III gives the comparison of the quantities e~ A ^ and B{x) obtained from 
least-square method applied to \og{h(d;x)/CiP(d)) vs d and analytical expressions 
for them. 

As a corroboration of the above conjectures we will obtain the formula for the 
maximal gap K(x) between two consecutive values of k < \fxj2 giving the prime 
4/c 2 + 1. Assuming, that the maximal gap K(x) appears only once we have the 
equation h(K(x),x) = 1 and putting the Hardy-Littlewood formula for ir q {x) in 



(70) and replacing P(d) by s we obtain: 



K(x) ~ ^ Q log(aO + log(2Cj) - 2 log(log(z)) ) (71) 



what for large x goes to the 



K(x)~—\og\x). (72) 



Because for all primes it is widely believed that G(x) = max( PniP „_ 1<a: )(p 

n Pn— l) ~ 

log 2 (x), we see that K{x) differs from G(x) just by a constant (but here only a frac- 
tion of k's are primes!). The largest gap between adjacent k giving 4/c 2 + 1 < 10 20 
was 290. The comparison of this formula with computer data is shown in the Figure 



6. From (72) we deduce the following 



Conjecture 4: 

q n+ i - q n = C( v /g^log 2 (g n )) (73) 

Let us recall that for all primes the Riemann Hypothesis gives p n +i — Pn = 0(pn +e ) 
for any e > 0, but in reality gaps between consecutive primes are smaller and 
the Cramer conjecture [10] states that p n +i — Pn = 0(\og 2 (p n )), see however [To] . 
Because our Conjecture 4 is obtained from the guessed formula for maximal gap 



between k's we expect (73) to be close to the optimal bound for g n+1 — q n . 
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7 Analog of the Chebyshev's bias 



7 Analog of the Chebyshev's bias 



For ordinary primes the Dirichlet's Theorem on the primes in arithmetical progres- 
sions asserts that the number 7r(x;4, 1) of primes < x giving 1 as the remainder 
when divided by 4 should be equal to the number 7r(x;4,3) of primes < x giv- 
ing 3 as the remainder when divided by 4, see e.g. [34] • But the direct inspec- 
tion shows that for small x there are more primes p = 3 (mod 4) than p = 1 
(mod 4): 7r(x;4, 3) > ir(x;4:, 1), what is called Chebyshev's bias [21]. For the first 
time l's takes the lead at P2946 — 26861 — up to this prime 3's win or there is a tie: 
ir(x; 4, 3) > n(x; 4, 1) for x < 26861. The next time n(x; 4, 3) < n(x; 4, 1) at 616841 
and in general there is preponderance of primes in the progression 4n + 3. The 
same phenomenon was observed in other arithmetical progressions [17] ■ However 
for the case of number of twins 7r(2; x) (primes pairs p, p + 2) and number of cousins 
7r(4;x) (primes pairs p, p + 4) which by the B conjecture of Hardy and Littlewood 
(35) should be the same ir(2; x) ~ 7r(4; x), there is no Chebyshev bias at least up to 
2 ?a 1.4 x 10 12 : sometimes twins and sometimes cousins take the lead, |41j . In this 
paper it was shown numerically that 7r(2;x) — 7r(4;x) behaves as the uncorrelated 
random walk; also the number of returns of this random walk to the origin (the 
number of such x that 7r(2; x) = 7r(4; x)) follows the usual square root law yfx [11. 
It is in agreement with last sentences of the paper [T7] suggesting that there is no 
Chebyshev bias for pairs of primes p,p + 2k in general. 



TABLE IV 



X 


7r q (x;3,2) 


7Tq(a;;3,2)/7r 9 (x;3, 1) 


W 2 (x)/ir q (x) 


Si(x) 


S 2 (x) 


8 (x) 


10 3 


8 


4.000000 


0.818182 


0.000000 


0.850210 


0.000000 


10 4 


12 


1.714286 


0.700000 


0.023927 


0.873320 


0.024337 


10 5 


32 


1.684211 


0.269231 


0.217321 


0.709371 


0.030297 


10 6 


71 


1.731707 


0.132743 


0.340751 


0.597523 


0.028523 


10 7 


215 


2.128713 


0.470032 


0.341211 


0.592423 


0.039598 


10 8 


560 


1.992883 


0.769596 


0.301877 


0.640202 


0.034797 


10 9 


1565 


1.924969 


0.322825 


0.363780 


0.582380 


0.033401 


10 10 


4410 


1.963491 


0.140454 


0.424147 


0.526904 


0.030546 


10 11 


12503 


1.978636 


0.049673 


0.476509 


0.479021 


0.027770 


10 12 


36069 


1.999279 


0.073571 


0.514925 


0.444198 


0.025570 


10 13 


104214 


2.009254 


0.671500 


0.476702 


0.485462 


0.025570 


10 14 


304482 


2.004754 


0.887649 


0.442651 


0.522215 


0.023710 


10 15 


893184 


1.999557 


0.610524 


0.442717 


0.524391 


0.022016 


10 16 


2636218 


2.000222 


0.709507 


0.433234 


0.535911 


0.020650 


10 17 


7817678 


1.999806 


0.474900 


0.442119 


0.528811 


0.019378 


10 18 


23265155 


1.999574 


0.169032 


0.471831 


0.500710 


0.018267 


10 19 


69497579 


1.999852 


0.141396 


0.491985 


0.481998 


0.017256 


10 20 


208240005 


2.000040 


0.365513 


0.499533 


0.475749 


0.016352 



All the primes of the form q = 4k 2 + 1 are necessarily disposed somewhere in 
the arithmetical progression 4/ + 1: q = 1 (mod 4). But the situation changes 
when we consider residues modulo 3. Let vr q (x; 3, 1) denote the number of primes 
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q = m 2 + 1 < x such that q = 1 (mod 3) and let 7r g (x; 3, 2) denote the number 
of primes q = m? + 1 < x such that q = 2 (mod 3). The direct inspection of all 
possibilities shows that the residue 2 should appear twice as often as the residue 1: 



It translates into the human 10-base system as the observation, that except for the 
first two cases 2 and 5, the last digit of the primes of the form m 2 + 1 can be only 
1 or 7, see ([I]). The direct inspection of all 10 possibilities under the assumption of 
normality of digits of k in the base 10 shows, that there are two times more ways of 
obtaining 7 than 1. 

The ratio of the number of those primes p = 4k 2 + 1 congruent to 2 modulo 3 
to those congruent to 1 modulo 3 can give some information about irregularities in 
the distribution of primes of the form 4k 2 + 1. The Table IV shows the values of 
the ratio n q (x; 3, 2)/ir q (x; 3, 1) for x = 10 3 , . . . 10 20 . As it is seen from this Table 
after the initial transient interval below 10 7 this ratio begins to oscillate around the 
predicted value 2. Initially ir q (x; 3, 2) > 2n q (x; 3, 1) and for the first time Ti q (x; 3, 2) < 
2n q (x; 3, 1) at qn = 4-42 2 + l = 7057, i.e. at the 17-th prime of the form q = m 2 + l. 
Recently A. Granville and G. Martin [17] have discussed several examples of "prime 
races". Primes of the form q = 4k 2 + 1 provide another example of such a race. 
Namely we will say that two wins at a given x if n q (x; 3, 2) > 2ir q (x; 3, 1) and let 
W2(x) denote the number of those primes q = 4k 2 + 1 < x that residue two wins over 
residue one. In the third column of the Table IV the ratio W2(x)/7r q (x) is shown. 
As it is seen from this sample of numbers there are large fluctuations of the ratio 
W2(x)/ir g (x). The Fig. 7 shows the plot of 7i q (x;3,2) — 2ir q (x;3, 1) up to x = 10 20 . 
This plot can be interpreted as a kind of one dimensional random walk: let y(x) 
denote the displacement of the walker at the "time" x, which plays the role of the 
time. If for a given q G Q we find that q = 1 (mod 3) the random walker performs 
step down of length 2 and if q = 2 (mod 3) the random walker performs step up 
of length 1 at the moment x = q. In other moments of time x the walker simply 
does not move. Thus we have y(x) = n q (x; 3, 2) — 2ir q (x; 3, 1). This plot resembles 
usual random walk, there were 21349 returns to the origin of this random walk up 
to 10 20 . There are large regions that y(x) > as well as y(x) < suggesting that 
there is no Chebyshev bias in the distribution of primes q n . In the Fig. 8 the plot of 
y(x)/x l l 4 = (vr q (x; 3, 2) — 2n q (x; 3, l))/^ 1 / 4 is shown. The amplitude of oscillations 
in this plot seems to be constant over the interval (10 3 , 10 20 ) and contained in the 
very short interval (— |, \). Dividing (ir q (x;3,2) — 2ii q (x;3, 1)) by x a results in 
amplitudes going to zero when a > \ and increasing when a < |. It is another 
argument in favor of the error term conjectured in the Sect. 2. There are probably 
logarithmic factors present, like for the usual Chebyshev bias, see Fig. 6 in [17] . but 
we are not able to separate it. The amplitude of oscillations of y(x)/x 1 ^ is very 
small, less than 0.5 and roughly half of the plot in Fig.8 is greater than zero and 
roughly half is below line zero. In [34] it was proposed to use the logarithmic density 
to measure the Chebyshev bias. Here we will define these densities for primes from 
Q as follows: 



7T g (x;3,2) w 2n q (x;3,l). 



(74) 



8 



i — 




2<n<x 
27r 9 (n;3,l)>7r 9 (n;3,2) 




n 



1 



(75) 
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7 Analog of the Chebyshev's bias 



lim 



lim 



log(z) 



log(x) 



E 



2<n<x 
27r,(n;3,l)<7r ? (n;3,2) 



?? 



E 



2<n<x 
27r 9 (n;3,l)=7rq(n;3,2) 



n 



(76) 



(77) 



We do not have at our disposal any formulas like those in [31] and we have to turn 
to the brute force numerical calculation of finite size approximations 5\(x), 82 (x) and 
Sq(x) given by expressions (75) — (77) without limit operation linL E _ KX) . The results 



are presented in the 5-th, 6-th and 7-th column of the Table IV and in Figure 9. Up 
to x = 2 31 the data for the Table IV and Figure 9 was obtained by direct summing of 
the harmonic sums, for x > 2 31 w 2.15 x 10 9 the incredible accurate approximation 



H, [ZD pp. 76-78] 



E 

k=n 



log ( m + 



log 



77 



+ 



(7? 



was used, thus the error of each summand was smaller than 10 -19 , and as there were 
(9(10 8 ) terms, in the worst case of adding up all roundoffs we expect the total error 
to be smaller than 10~ 10 . To calculate the harmonic series up to x = 10 20 directly 
by adding all numbers \jn would take a few thousands years of CPU time, but it is 
in general impossible using standard programming languages as the loop can have 
only integer counter and on 64 bit processors the largest integer is 2 63 9.22 x 10 18 . 
For all primes Chebyshev conjectured that 



lim £(-1) 



-00 



(79) 



p>2 



It was proved by Hardy and Littlewood [18] and Landau [27] , [26] that ( 79 ) is true 
if and only if the L-function 



L(s,xi) = 



n=0 



(2n + 1) { 



does not have nontrivial zeros outside the critical line 3?(s) 
the analogous function for the case of primes from Q: 



JO) 



Here we formulate 



where 



F(x) 



2 

-1 



E 



c q e 



-q/x 



when q mod 3 
when q mod 3 



(82) 



We made the plot of F{x) and in contrast to (79) this function seems to not possess 
a limit when x — > 00. Because values of F(x) sometimes are close to zero calculating 



the sum (81) on the computer we have stopped summation at such q' that 



eye 



-q'/x 



c„e 



■2 ^1 



-q/x 



< 10 



-8 



(83) 
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thus we have used relative error. Using such a condition (83) for terminating the 



sum in (81 ) is necessary as F(x) sometimes crosses zero and using absolute error can 
be misleading at small values of F(x). The last 43 points in the Fig. 10 were not 
fulfilling the requirement (83) as all primes q = m 2 + 1 generated were exhausted, 



but as it is seen in the plot of F(x) in Fig. 10 values of the function F(x) are well 
above 1000 for x > 10 17 . 

The numbers presented in the Table IV and plots in the Figures 8, 9 and 10 
allow us to formulate the 

Conjecture 5: There is no Chebyshev bias for primes of the form m 2 + 1: 
$l = 8 2 = §, 5 = 

It is the last conjecture formulated in this paper. 
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Fig.l The plot of the error term for the conjecture Q up to x = 10 20 on the 



tt„(x) ' ' ,hi 



in r 

2^1 J2 



was 



double logarithmic axes. The plot of the difference , ,v ,, ,,, ^ , 

made from the 101566 points. Up to 10 12 values of |A 9 (x)| at all 54110 primes of 
the form m 2 + 1 < 10 12 are included. For larger values of x > 10 12 the following 
decimation procedure was used: changes of local maxima larger than 1, all local 
minima and all sign changes were recorded with the provision that values 
corresponding to the sign changes of the |A ? (x)| and smaller than 10~ 3 were 
artificially set to 10~ 3 . Additionally the plot contains 18430 points recorded at the 

geometrical progression x = 10 12 x (1.001)™. In blue the power-like fit 
0.5 1 877 x x 2139 to u(x) obtained by the least-square method is shown and in red 
the possible choice for the bound of oj(x) from above is plotted. 




Fig.2 The plot of the oscillations of the (ir q (x) — \C q ^ d 1 " gM )/x 1//4 . It consist of 
106027 points and the procedure for recording values of this difference was similar 

to the previous one. 
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Fig. 3 The plot of the number of sign changes of the difference 
7r ? (x) — \C ^Aog{u) U P to x = 10 20 . There are 20456 data points and power-like 
fit was performed with respect to all points. In red the power fit obtained by the 

least-square method is shown. 
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Fig.4 The plot of the function ir q (d; x) for x = 10 10 , 10 12 , . . . , 10 20 and d < 300. By 
red boxes the values of n q ( y d;x)/P(d) are shown; for clarity only every 4-th value of 
d is shown — descarding small fluctuations they fit the green lines representing the 



integrals C\ f r 



dt 



5 v^logW 
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Fig.5 The plot of the histogram h(d; x) for x = 10 10 , 10 12 , . . . , 10 20 . . By red boxes 
the values of h(x, d)/P(d) are shown; for clarity only every 4-th red box is shown. 
They are supposed to lie along the green lines representing the plots of the 





Fig.6 The plot of the function K(x) 
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Fig. 7 The plot of the difference (7r 9 (x; 3, 2) — 27r q (x; 3, 1)). There are 102714 points: 

up to 10 12 all primes q n are plotted, for x > 10 12 the following decimations 
procedure was used: all local maxima, minima and sign changes of this difference 

were recorded plus 18430 values at the geometrical progression 10 12 x (1.001) n . 
Because the amplitude grows with x we have presented the interval (10 10 , 10 20 ) in 
the inset with smaller range on the y-axis. 




Fig.8 The plot of the ratio y{x)/x 1 ^ = (7r 9 (o;;3, 2) — 27i q (x; 3, l))/x 1 ^. This plot is 
made of 102714 points, exactly as previous plot. 
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Fig. 9 The plots of the logarithmic densities Si (x) , 82 {%) and 5q(x) defined in the 
text. Each plot consists of 72542 points: up to 10 12 all all primes q n are plotted, 
for x > 10 12 the values of <5(x)'s were recorded at the progression 10 12 x (1.001)™. 
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Fig. 10 The plot of the function F(x) defined by equation (79). There are 2021 

points on these figures: the arguments x were chosen as the geometrical 
progression 100.0 x 1.02) n , n — 0,1, ... , 2020 and the largest value of x is here 
2.35693149 . . . x 10 19 , for which the last term in the sum (81) was roughly 

e -10 2 °/2.35693149xl0^ = Q.01436723 .... 



